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ABSTRACT 

Simple stellar population (SSP) synthesis models are useful tools for studying the nature of 
unresolved star clusters in external galaxies. However, the plethora of currently available SSP 
models gives rise to significant and poorly documented systematic differences. Here we con- 
sider the outputs of the commonly used Bruzual & Chariot and GALEV models, as well as 
a recently updated SSP model suite which attempts to include the contributions of binary 
merger products in the form of blue straggler stars (BS-SSP). We rederive the ages, metallic - 
ities, extinction values and masses of 445 previously observed globular-like clusters in M3 1 
based on \ 2 minimisation of their spectral energy distributions with respect to these three 
different SSP models and adopting a Chabrier-like stellar initial mass function. A compari- 
son between our new results and previous estimates of the same parameters shows that the 
Bruzual & Chariot models yield the youngest ages and lowest masses, while adoption of the 
BS-SSP models results in the oldest ages and highest mass estimates. Similarly, the GALEV 
SSP models produce the lowest metallicities, with the highest values resulting from the BS- 
SSP model suite. These trends are caused by intrinsic differences associated with the models, 
and are not significantly affected by the well-known age-metallicity degeneracy. Finally, we 
note that the mass function of the massive M3 1 star clusters is similar to that of the Milky 
Way's globular clusters, which implies that the two star cluster systems likely formed under 
similar environmental conditions. 
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1 INTRODUCTION 

Star clusters represent an important stellar population compo- 
nent of galaxies. Their age distributions trace the main evolution- 
ary events associated with the formation and evolution of their 
host galaxies. However, most clusters in external galaxies can- 
not be resolved into individual stars, even with the superb reso- 
lution of the Hubble Space Telescope (HST). Consequently, stel- 
lar population synthesis has become an important and power- 
ful tool to interpret the nature of extragalactic star clusters, in- 
cluding their ages, metallicities, reddening values and masses. 
In the past few decades, many different simple stellar popu- 
lation (SSP) synthesis models have been constructed and ap- 
plied to study extragalactic s tar clusters based on phot o metry in 
multip l e passbands, inc luding Bruzual & Chariot ( 1993); Worthevl 
dl994h : iLeitherer & Heckmanl d 19951) : iFioc & Rocca-Volmerangd 
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); Bruzual & Chariot 


(2003 


);lYietal. 1 


2003); 


Maraston (2005); 


Conrov et al.l d2009); 


XinetalJ < 201 lb. 


as well as the GALEV models (see, e.g., 


Schulzetal. 2002 


L 2003; Anders & Fritze- v. Alvenslebenl 120031; 



cently updated BS-SSP models jXin et alj201 If) considered the ef- 
fects of blue straggler stars (BSs), which are likely the products 
of binary mergers and which seem a common occurrence in star 
clusters ( see, e .g., lAhumada & Lapassedfl995l ; IPiotto et al.l [20021; 
IXin et alj|201 1). Since BSs are often luminous and can render the 
integrated luminosities and colours of their host clusters signifi- 
cantly brighter and bluer, we set out to explore the differences 
between results derived from application of the BS-SSP models 
and those derived fro m other models. Here we choose two com- 
monly used models - Bruz ual & Chariot! d2003l henceforth BC03) 
and GALEV - for our comparisons. 

Star clusters were long thought of as members of two distinct 
types, i.e., open clusters (OCs) and globular clusters (GCs). OCs 
are generally young, not very massive, faint, diffuse, and usually 
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located in galactic discs, quite contrary to the nature of GCs, which 
are mostly old, massive, luminous, centrally concentrated, and usu- 
ally located in the haloes of their host galaxies. However, this sim- 
plistic picture has been changing since the discovery of young 
massive star clu sters (YMCs) in man y galaxies, in cluding in the 
Milky Way (e.g..lAscenso et alj2007allbl), M3 1 (e.g..lBarmbv et all 



200' 



9; 



2009 I- ICaldwell et alj 

20 1 ;lHodge et alj2010h. M82 (e.g..lO'Connell & Mangan 



Ma et al 

T~ 



2009; Perinaetal 



2009 



1978 



de Grijs et al. 2003a; Smith et al. 2007; McCrady 2009), the An- 



tennae galaxies (e g..lWhitmore et al.ll2005l) and NGC 11 40 (e.g., 



iHunter et alJI 19941 ; Ide Grijs et al.ll2004 iMoll et alj|2008l) as well 
as NGC 3603 and R136 dMafz Apellanizll2009l) . YMC properties 
span those of both OCs and GCs, with typical masses (> 10 4 
M0) greater than those of (most) OCs and young ages (< 1 Gyr), 
i.e., quite different from present-day GCs, so that they are often 
considered candidate proto-GCs. The new category of YMCs ren- 
ders cluster classification blurred and difficult. In this paper, we use 
the term 'globular-like' cluster to distinguish massive (YMCs and 
GCs) from lower-mass clusters (OCs). Since OCs are usually faint 
and located in galactic discs, this makes them difficult to study, 
which is why we focus on globular-like clusters here. 

Located at a distance of ~780 kpc dstanek & Garnavichl 
1 19981 ; lMacrill200ll ; iMcConnachie et al.ll2005f) . M31 is the nearest 
large spiral galaxy in our Local Group of galaxies. Therefore, it 
represents an ideal laboratory for studies of statistically signifi- 
cant numbers of globular-like clusters in external galaxies. In ad- 
dition, since the Hubble type and mass of M31 are similar to those 
of our Galaxy, studying M31's massive cluster system is also im- 
portant for our understanding of the Galactic GC system. Based 
on HST Wide Field an d Planetary Camera-2 (WFPC2) images, 
iKrienke & Hodgd 120071) suggested that there may be ~ 80,000 star 
clusters in the M31 disc. Most of these disc clus ters are faint OCs. 
The n umber of GCs in M31 is much smaller. iBarmbv & Huchra 



2001) estimated their total number at 460 ± 70, while Perinaetal 
2010h arrived at ~530, with an additional ~100 YMCs. To limit 



the scope of this paper, we only focus on the globular-like clus- 
ters, including GCs and YMCs, because they are luminous and 
relatively easy to observe at the distance of M31. Studies aimed 
at identification, classification and analysis of the population of 
M3 1 globul ar-like clusters have been undertaken since the pioneer- 
ing work oflHubbld dl932l) ( s ee, e.g ., Vetesnik 1962; Sarge nTet al.l 
[19801 " 



19771 IBattistini et al.l [198 
Barmbv et alfcOOOh 



Il987l. Il993l ; ICrampton et alj Il985l : 



Application of the \ -minimisation technique to spectral- 
energy-distribution (SED) fitting is a commonly used method for 
estimating ages, metallicities, red dening values and masses of ex- 
tragala ct ic star clusters (see, e.g.. I jiang et a l. 2003; d e Grijs et alj 
2003allbl 120051; Ide Grijs & AnderslboOol iFan et alj|2006t iMa et ail 
20071 f2009; Wang et al. 201Ch . The M31 massive star cluster sys- 



tem has also been studied extensively in this manner. For instance , 
I Jiang etalJd2003h studied 172 GCs from the lBattistini et alldl987l) 
sample which were located in the central ~ 1 deg 2 region of the 
image and observed by the Beijing-Arizona-Taipei-Connecticut 
(BATC; see IFan et alj 19961) survey. The observations of all of these 
clusters had relatively high signal-to-noise ratios (S/N) - and thus 
good photometry - in the 13 BATC intermediate-band filters cov- 
ering the range from ~ 4, 000 to ~ 10, 000 A. By comparison of 
their SEDs with the SSP mo dels of Bruzual & Chario t (1996; un- 
publis hed, updated version o flBruzual & Charlodl993l) . |jiang et"al] 
d2003l) concluded that nearly all their sample GCs were older than 
1 Gyr. (Note that their sample GCs were mainly selected from 
the M31 bulge and inner disc, while they only used three differ- 



ent metallicities for their comparisons.) Fan et al. I d2006h estimated 
the ages of 91 GCs in M31 from the ljiang et alj [2003) sample by 
matching BATC intermediate-band and Two-Micron All-Sky Sur- 
vey (2MASS) JHK SEDs with BC03 SSP models. They identified 
a young, ~3 Gyr-old population and an intermediate-age popula- 
tion of ~8 Gyr, in addition to the well-known old (>10 Gyr) GC 
population. 

It has been unequivocally established that near-infrared (NIR) 
photometry can partially break the a ge-metailicity degen e racy and, 
hence , improve age estimates (e.g., Anders et al.ll2004bT) . IMa et al.1 
d2007h determined the age of an old M31 globular-like clus- 
ter (S3 12) based on GALEX near-ultraviolet, optical broad-band, 
BATC and 2MASS JHK photometry. They concluded that this 
cluster has a mass of (9.8 ± 1.85) x 10 s M Q and an age of 9.5 
Gyr. Subsequently, IMa et alj d2009h used the GALEV SSP mod- 
els applied to BATC, 2MASS JHK and GALEX data to derive 
the ages of 35 GCs in th e central M31 field that were not in- 
cluded inl Jiang et ail ^ 2003). These clusters were also selected from 
IBattistini et alj ( fl987l) and located in the galaxy's central ~ 1 deg 2 , 
but most of their observations were characterised by relatively poor 
S/Ns. (Note that the GALEV SSP models provide spectral templates 
for ages up to 16 Gyr and metallicities Z > 0.0004, while the 
BC03 models include ages up to 20 Gyr, so that minor differences 
may result for the oldest and lowest-metallicity (Z < 0.0004) clus- 
ters, depending on the SSP model suite adopted.) IMa et ai] d2009h 
found that most of their sample GCs cover ed the age range f rom 1 
to 6 Gyr, with a peak at ~ 3 Gyr. Recently. IWang et al. (2010) per- 
formed photometry of another 104 M3 1 globular-like clusters using 
BATC multicolour observations of the central ~ 6 deg 2 M3 1 field, 
which covers a large part of the galactic disc, where many young 
star clusters can be found. They estimated the clusters' ages by fit- 
ting their SEDs with GALEV SSP models, revealing the presence of 
young, intermediate-age and old cluster populations in M3 1 . All of 
these studies were based on SED fitting of multicolour photometry 
and comparison with SSP models. 

Despite a significant historical body of work and the large 
number and diversity of currently available SSP models, sys- 
tematic differences among model outputs persist. For instance, 
IPeacock et al.l d201ll) recently compared several SSP models with 
photometry obtained through the Sloan Digital Sky Survey filters of 
M31 globular-like clusters. They identified a signific ant offset be- 
tween the models and the observed (g — r) colours. IConrov et al.l 
(2009) investigated the propagation of uncertainties in stellar popu- 
lation synthesis modelling, and specifically their impact on our un- 
derstanding of the observational manifestations of stellar evolution, 
the stellar initial mass function (IMF) and the luminosity evolu- 
tion of stars and stellar populations. They found that stellar masses 
could be affected by errors of ~ 0.3 dex in nearby galaxies, while 
for more distant bri ght red galaxies, the unc ertainties in mass are at 
the ~ 0.6 dex level. lconrov & Gunnl d2010h compared and assessed 
the differences in the luminosities and colours of their stellar pop- 
ulation synthesis models (FSPS) with r espect to those associated 
with both the BC03 and Ma rastonl d~2005l) (M05) model suites. They 
found that the FSPS and BC03 models perform well for star clus- 
ters in the Magellanic Clouds, while the M05 models seem to be 
too red and result in incorrect age dependences. However, all three 
models provide poor fits to the NIR colours of star clusters in both 
the Galaxy and M31. They also found that the FSPS models can 
fit ultraviolet photometry successf ully, while the BC03 and M05 
models fail in doing so. Recently, iMaraston & Strombackl d201lh 
compared the colour indices and spectral profiles of their models 
based on application of different empirical stellar spectral libraries, 
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including the observational Pickles, ELODIE, STELIB and MILES 
libraries, the theoretic al MARCS stellar library, a nd the models of 
IVazdekis et ail ( bold) and lConrov & Gunnl fcOlCh . They noted that 
SSP models based on empirical libraries have lower fluxes in the 
V band than theoretically predicted by t he BaSeL spectral library 
dLeieune. Cuisinier & BuseJl997l[T998l) , which in turn is based on 
the lKuruczld 1992b library. This could lead to bluer (B — V) colours 
(by up to 0.05 mag) in their models. In this paper, we aim at explor- 
ing the causes and extents of systematic effects resulting from ap- 
plication of specific sets of SSP models to broad-band photometric 
observations. This paper builds on our previous work in characteris- 
ing uncertainties o wing to application of SSP models to broad-b and 
SEDs, in particularlAnders et al.ld2004bl) . Ide" Griis et al l d2005h and 
Ide Griis & Anders! J2006I) . Although a small fraction (< 5%) of our 
sample clusters have masses < a few x 10 4 M©, we specifically do 
not consider the effects of stochastic sampling of the stellar IMF. 
Stochasticity will increase the uncertainties associated with the re- 
sulting model parameters, particularly so for the lower-mass clus- 
ters, i.e., < a few x 10 4 Mq (e.g., Cervino, Luridiana & Castander 
2000; Cervino et al. 2002; Cervino & Luridiana 2004, 2006; Mai'z 
Apellaniz 2009; Popescu & Hanson 2010; Fouesneau & Lancon 
2010; Fouesneau et al. 2012; Popescu, Hanson & Elmegreen 2012), 
but a detailed study of these effects in the context of our current 
work is beyond the scope of this paper. 

Here we derive the ages, metallicities, reddening values and 
masses of the 445 confirmed globular-like clusters and cluster can- 
didates by comparing the observed SE Ds with GALEV models as 
well as the latest BS-corrected models dXin et al.ll20~ lh since re- 
cent studies show t hat nearly all star clusters contain so-called blue 
straggler stars ( e.g., lAhumada & Lapassetlll995UPiotto et alj2002l ; 
IXin et al]|201lh . which could artificially render the integrated star 
cluster colours much bluer than their real, intrinsic colours. It is 
interesting to compare the results from different models, which 
helps us to better understand the practical impact of the differ- 
ences between the models. In addition, we compare our BS-SSP 
and GALEV-based results with those derived using the BC03 model 
suite. This paper is organized as follows. Sect.|2]describes our sam- 
ple selection and UBVRIJHK photometry. In Sect. [3] we de- 
scribe how we fit the SEDs using the different SSP models. Sect. 
|4]presents the resulting ages, metallicities, masses and their distri- 
butions as well as comparisons of the estimates derived based on 
adoption of different SSP models. Finally, we summarise and con- 
clude the paper in Sect. [5] 



2 SAMPLE 



iFan et al. I d2010h provided the basic globular-like cluster sample 
and photometry for our work. They performed UBVRI photom- 
etry based on archival images of the Local Group Galaxy Sur- 
vey, which covers a region of 2.2 deg 2 along M31's major axis. 
The fields were observed from August 2000 to September 2002 
with the KPNO M ayall 4 m telescope (for observati onal details see 
iMassev et al.ll2006h . The coordinates specified by ICaldwell et al.l 
d2009h were used for cluster identification and aperture photometry 
based on a set of variable aperture radii ranging from 1.03 to 9.06 
arcsec. A comparison of this new photometry with that published 
in the Revised Bo logna Catalogue (RBC) v.4 jGalleti et a"l]|2004 
2006 1 l200l 120091), as wel l as with other data sets dBarmbv et al.l 
2000; Peacock etal.ll2010h . showed minor to negligible system- 



their sample. lFan et al.l d2010h also included a small fraction of the 

UBVRI photometry from the RBC v.4. Thus, the UBVRI cluster 

I 1 J — ~~ 71, 

photometry adopted in this paper is homogeneous. Fan et al. ( 2010) 
also adopted the JHK photometry of lGalleti et al.l a2009h . because 
NIR photome try is important to ( partially) break degeneracies in 
SED fits (e.g., lAnders etail2004bl) . Finally, the authors obtained a 
sample of 445 confirmed globular-like clusters and cluster candi- 
dates with photometry in at least six of the U BVRIJHK bands. 
They derived the ages, metallicities, masses and reddening values 
based on the BC03 model suite, using Pado va isochrones an d a 
Chabrier IMF. Therefore, we use the sample of lFanetal.ld2010l) as 
the basis for our tests and comparisons. 



3 METHOD 

In this section, we describe the methods and processes used for 
our redeterminations of the cluster ages, metallicities, extinction 
values and masses based on SED fitting. T he gen eral fitting method 
used here is the same as that in lFan et aU d2010f) . although we use 
different, updated SSP models. 



3.1 SSP models used 

Since we h ave access to pho tometry in optical and NIR broad-band 
filters from lFan et al. Id20l3), we can obtain the clusters' fundamen- 
tal integrated physical parameters, such as their ages, metallicities 
and masses, by means of SED fitting. However, a large body of 
evidence suggests that a strong age-meta lficity degeneracy domi- 
nates if only optical pho t ometry is us ed (Worthevl ll994t [Arimoto 
ll996l ; lKavirai et alj2007l) . lAnders et al] d2004bl) concluded that the 
degeneracy can be partially broken by adding NIR photometry to 
the optical colours, with the ef ficacy of this process dep ending on 
the ag e of the stellar population. Ide Griis et al.l J2005h and lWu et al.l 
(2005) also showed that use of NIR colours can greatly contribute 
to break the age-metallicity and age-extinction degeneracies. It has 
also been shown that [/-band data is important to obtain accurate 
age, metallicity and ext i nction es timates based on SED fits (see, 
e.g.jAnders et al. 2003a, 2004a bl;lde~Griis et al.l2005l;lGieles etall 
l2005l ; lBastianeta l. 2005a d lKotulla et alj|2009l). Thus, in our fits , 
we used the full U BVRI photometric data set of lFan et all d2010h , 
complemented with JHK photometry from the RBC v.4, to disen- 
tangle the degeneracies and obtain more reliable results. 

We used a \ 2 -minimisation technique for our age, metallic- 
ity, mass and extinction estimates, comparing the observed, inte- 
grated SEDs with theoretical SSP models (see, e .g. jFan et al.l2006l . 
|2010| ; iMa et alj|2007l , 120091 ; I Wang et al.ll201(il) . In this paper, we 
take advantage of the availability of a new set of SSP models tha t 
attempt to correct for the contributions of BSs dXin et al.l |20 lib . 
The BS-SSP model suite is based on the Padova 1994 isochrones 
and the BaSeL-Kurucz stellar spectral library. We decided to use 
these BS-SSP models because the (central) cluster colours are of- 
ten affected by the presence of BSs, thus causing deviations in 
(predominantly) the U and B bands with respect to SSP models 



atic differences, with maximum offsets of 0.1 to 0.2<r, where a 
refers to the r.m.s. scatter in the photometric differences. To enlarge 



based solely on single-star evolution dXin & Deng|l2005l ; IXin et ail 
120071 ; ICenarro et al]|2008l) . Following lxin et al.l d201 lb. w e adop t 
a t wo-part power- law IMF, similar to both the IChabrierl d2003l) 
and lKroupal d200ll) formalisms. The BS-SSP synthesis models in- 
clude five initial metallicities {Z = 0.0004, 0.004, 0.008, 0.02 = 
Z Q , and 0.05, (corresponding to [Fe/H] = -1.7,-0.7,-1.4,0 
and +0.4 dex), and 24 equally logarithmically-spaced time steps 
from 100 Myr to 20 Gyr with a wavelength coverage from 91 A 
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to 160 11m. To avoid discontinu i ties, followi ng Fan et al. 1 20061). 
iMa et alj l2007l) . lMa et alj l2009h . lWang et alj feOld) and lFan et all 
d2010h . we linearly interpolated the input metallicities to create a 
higher-resolution spectral grid containing 100 metallicities (from 
[Fe/H] = -1.7 to +0.4) and 231 ages from 0.1 Gyr to 20 Gyr, 
with equally logarithmicall y-spaced intervals betwe en the newly 
created templates. Although lFravn &"G ilmore (2002) discussed the 
merits and problems associated with linear interpolation between 
two points in logarithmic space, wh ich may especially be an issue 
for sharp 'hooks' in our isochro nes, ICharlot & Bruzuall(ll99ll) and 
Bruzual &Charlotl ( fT993l l 20031) suggested that interpolation of the 
model ages and metallicities is feasible and can be done reasonably 
well. 

The GAL AXE V models (BC03) follow the spectral and pho- 
tometric evolution of SSPs for a wide range of stellar metallici- 
ties. The model suite includes 26 SSP models at b oth high and low 
resolution, 13 of which were computed using the IChabrierl ( 2003 ) 
IMF assuming lower and upper stellar mass cutoffs of triL = 0.1 
and mu — 100 Mp , respectivel y. The other 13 models were com- 
puted using the lSalpeterUl955l) IMF with the same mass cutoffs. 
This model suite provides SSP models based on both the Padova 
1994 and 2000 evolutionary tracks. However, the authors point 
out that the Padova 2000 models tend to produce worse agree- 
ment with observed galaxy colours. These SSP models contain 221 
ages, in unequally-spaced time steps from 1 x 10 5 yr to 20 Gyr. 
The evolving spectra include the contribution of the stellar com- 
ponent at wavelengths from 91 A to 160 /im. In this paper, we 
adopt the high-resolution SSP models computed using the Padova 
1994 evolutionary tracks and a Chabrier (2003) IMF. The spec- 
tral libraries used include the theo retical Ba S eL a nd observational 
STELIB and Pickles collections. iFan et al linearly inter- 

polated the metallicities provided by the Padova 1994 model grid 
([Fe/H] = -2.2490, -1.6464, -0.6392, -0.3300, +0.0932 and 
+0.5595 dex) using 100 equal steps and fitted the SEDs of the 
globular-like clusters based on the high-resolution BC03 models. 
In this paper, we do not redo the fits based on th e BC03 models, 
but we simply adopt the results of lFan et all J2010h . 

For the GALEV models, we used the Padova 1994 theoret- 
ical isochrones with a Kroupa IMF and the BaSeL library of 
stellar spectra. Indeed, BC03 showed that the spectral proper- 
ties based on a Kroupa IMF are very similar to those using the 
Chabrier IMF. The models include five initial metallicities, Z = 
0.0004,0.004,0.008,0.02 and 0.05, corresponding t o [Fe/H] = 
-1.7,-0.7,-0.4,0 and +0.4 dAnders et alj |2003bl) with 4000 
equally-spaced time steps from 4 Myr to 16 Gyr (in steps of 4 Myr). 
In fact, to approximate a comparable resolution in age, we reduced 
the number of time steps to 200: we kept the model's original time 
steps for ages < 1 Gyr and enlarged the time interval to 0. 1 Gyr for 
older ages. The wavelength coverage also spans the range from 91 
A to 160 /im. Similarly, to achieve higher metallicity resolution and 
avoid discontinuities, we enlarged the spectral grid - which con- 
tained 105 metallicities (from [Fe/H] = —1.7 to 0.4) - by inter- 
polating in logarithmic space, using equally logarithmically-spaced 
intervals between th e newly created templates, as recommended by 
lKotullaetalJj2009l) . 

Using Fig. Q] we investigate the UBVRIJHK colour be- 
haviour as a function of age as predicted by the different mod- 
els, i.e., the BS-SSP, BC03 and GALEV models. We adopted so- 
lar metallicity and all colours are in the Vega photometric sys- 
tem. The BS-corrected models yield offsets to the blue for all 
colours: for ages > 2 Gyr, we find colour offsets (U — E) ~ 0.2, 
(B-V) ~ 0.1, (V-R) ~ 0.1 and (V-/) ~ 0.1 mag with respect 



to the other models, while for ages < 2 Gyr, the offsets for these 
same colours are smaller or even negligible. However, the GALEV 
models exhibit much redder colours in ( V — K ) and (J — H) for 
ages < 2 Gyr, which may affect our age estimates of young stellar 
populations. The GALEV and BC03 models predict similar colours 
in (U - B), (B - V), (V - R) and (V - I) while the two model 
sets show significant offsets in (V — K) and (J — H) colours, 
particularly for young ages. 

Figure |2]is similar to Fig.[T]but shows the colours as a func- 
tion of metallicity for an age of 12 Gyr, where any differences are 
most prominent. In all panels, we find that the BS-SSP models pre- 
dict the bluest colours, while the GALEV models predict the reddest 
colours. Note also that for (V - R), (V - I) and (J - H), the 
GALEV and BC03 models predict similar colours (within reason- 
able observational uncertainties) as a function of metallicity. 

In Fig. [3] we show the mass-to-light ratios (M/Ls) in the V 
band as a function of age for solar metallicity for the three models. 
The GALEV models predict the highest M/Ls, while the BS-SSP 
models predict the lowest for all ages. In addition, the differences 
becomes larger for older ages. 

These intrinsic colour and M/L differences will lead to differ- 
ences in the results predicted by the different models. In particular, 
the BS-SSP models are expected to predict intrinsically older ages, 
higher metallicities and lower masses for the same age, regardless 
of the reddening affecting the cluster photometry. Fig. [3]shows that 
one should be aware that if the ages predicted by BS-SSP models 
are not sufficiently old, the M/Ly values derived from the BS-SSP 
models could be lower than those predicted by other models. 

3.2 Fitting procedure 

The BS-SSP and GALEV spectra can easily be convolved to mag- 
nitudes in the AB system using the filter-response functions in 
the UBVRIJHK bands. The apparent magnitudes of the BS- 
SSP/GALEV synthesis models in the AB system are given by 



tiab(£) = -2.5 log 



dA A R(\) 



48.60, 



(1) 



where -R(A) is filter-response function and Fa (A, t) the flux, which 
is a function of wavelength (A) and evolutionary time (t). Ai and 
A2 are the lower and upper wavelength cutoffs of the respective 
filter (see BC03). We converted all observed integrate d magnitudes 
(UBVRIJHK) to the AB system using fhe lKuruczl ll992f) SEDs. 
Since a comparison o f the exti n ction values of the M31 
lobul ar-like clusters from IFan et aL I d2008l) and iBarmbv et alj 
2000) shows that the offset is —0.01 ± 0.10, we continued on 
the assumption that there is no systematic offset between the two 
data sets. Therefore, for those ( 310 of 445) clust ers in our sample 
that ha ve reddening values from lFan et alj d2008l) or lBarmbv et alj 
( 2000 ), the magnitu des were corrected for reddening assuming a 
ICardelli etal.l ll989) Galactic extinction law, so that their ages and 
metallicities can be determined by comparison of the interpolated 
high-resolution BC03 SSP synthesis models with the SEDs from 
our photometry and leaving the metallicity [Fe/H] as a free param- 
eter, i.e., 



Xmin(t, [Fe/H]) = min 



8 

E 



(2) 



where m™ od (t, [Fe/H]) is the integrated magnitude in the i th fil- 
ter of a theoretical SSP at age t and for metallicity [Fe/H], m A bs 



Systematic differences in SSP models: M31 star clusters 5 




represents the observed, integrated magnitude in the same filter, 
m x . = UBVRIJHK, and 



2 

°obs 



+ On 



(3) 



where a hs,i is the observational uncertainty. Since the RBC 
does not include any magnitude uncertai nties, we applied the 
rough estimates from iGalleti et aD d2004h . i.e., 0.05 and 0.08 
mag uncertainties for the BVRI and U bands, respectively. As 
for th e NIR JHK m agnitudes, the uncertainties are estimated 
as in Fan et al ] d2006h by applying the relations in Fig. 2 of 
ICarpenter. Hillenbrand & Skrutskiel d200 ll) . which shows the ob- 
served uncertainties as a function of magnitude for b right stars in 
the 2MASS JHK bands. In addition. lFan et alj j2006) showed that 
the adopted uncertainty does not affect the quality of the SED fits. 
c"mod,i represents the uncer tainty associated wi th the model itself, 
for the i th filter F ollowinglde Griis et al.ld2005l) (the i r Sect 3.2.4) 
IWu et alj d2005b iFan et all d200d l201dh . iMa et alj d200l 120091) 
and lWang et al.l d2010l) , we adopt a mo d,i = 0.05. 

For the 135 of 445 sample clusters without extinction values 
from the literature, we constrained the ages, metallicities and ex- 
tinction values while keeping [Fe/H] as a free parameter, using 



xl in [t,We/R},E(B-V)}^mm 



E 



rn obs 



m mod " 



.(4) 



We varied the reddening between E(B — V) = 0.0 and 2.0 mag 
in steps of 0.02 mag. We then obtained the values for the extinction 
coefficient, R\, by i nterpolating the interstellar extinction curve of 
ICardelli et aL1dl989l) . We thus fitted the extinction-corrected SEDs, 
for which the model with the minimum \ 2 returned the best-fitting 

E{B - V) values. 

We used the same fitting method as in IFan et al] d2010h . If the 
initial age estimate of a given star cluster is older than the Wilkin- 
son Microwave Anisotropy Probe (W MAP) age of the Universe, 
13.7 ± 0.2 Gyr dSpergel et aT]|2003l) . we adopt an age of 12 Gyr 
as initial guess and iterate until the fitting routine reaches a lo- 
cal minimum. Although the estimated age could be older than 12 
Gyr, this approach e nsures tha t the e stimated age will not exceed 
the WMAP age: see IFan et afl fcOlCh for details. It is well-known 
that SSP SE Ds are not sens itive to changes in age for ages > 10 
Gyr (see e.g. jMa et al]|2007l and Fig. [TJ. Therefore, although the 
upper age limit in the BS-SSP and GALEV models is 20 and 16 
Gyr, respectively, the ages of the clusters determined here do not 
exceed the WMAP age (but see below for a discussion of the ef- 
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Figure 2. As Fig. 1, but for the predicted colours as a function of metallicity for 12 Gyr-old SSP models. 
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Figure 3. As Fig. 1, but for the predicted mass-to-light ratios in the V band as a function of age for solar metallicity. 
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fects of adopting an upper age limit in our fits). To estimate the 
la uncertainty of a given parameter, we first fix all other param- 
eters to their best-fitting minimum \ 2 values. We then determine 
Ax 2 = X 2 ~ Xmin < !• The confidence interval of this param- 
eter is then given by the la width of the normal distribution of 
Ax 2 , which thus yields the la uncertainty. We determine a clus- 
ter's mass based on the model's M/L for its specific age, metal- 
licity and de-reddened lumino sity. This is similar to the approach 
adopted and advocated by, e.g jAnders et al.ld2004bl) . Since 310 of 
445 clusters in our sample have robustly determined reddening val- 
ues available from literature sources, only 135 of the 445 clusters 
may suffer more significantly from degeneracies related to the un- 
certain amount of extinction affecting the clusters. However, we re- 
iterate that our NIR photometry is expected to enable us to partially 
break the age-metallicity degeneracy (but see Sect. 14. 1) . Therefore 
the age-metallicity-extinction degeneracy will most likely not be 
serious when we set the maximum model age to the WMAP age. 
In fact, for the BS-SSP models, 120 of our 445 clusters are older 
than the WMAP age, while for the BC03 and GALEV models the 
equivalent numbers are 55 and 99, respectively. 



4 FIT RESULTS AND COMPARISONS 

The fitting procedure outlined above allows us to obtain the ba- 
sic parameters of our sample clusters, i.e., their ages, metallicities, 
masses and reddening values. Therefore, we can now compare the 
fit results derived based on adoption of different models and differ- 
ent IMFs. 

4.1 Comparisons 

To investigate the differences caused by adopting free versus fixed 
reddenin g, we co llected data for 310 clusters with reddening val- 
ues from lFan et al.l ( 120081) and lBarmbv et al.l ( 120001) . and performed 
free- and fixed-reddening fits. We also wanted to check whether the 
free-reddening method could lead to degeneracies in our simulta- 
neous age, metallicity and mass determinations. 

Figure [4] shows comparisons of the extinction values, ages, 
metallicities and masses based on adoption of the BS-SSP models 
with reddening values from the literature, and those resulting from 
leaving reddening as a free parameter. Approximately half of the 
newly derived reddening values are systematically larger than the 
corresponding literature values. Similarly, half of the redetermined 
ages are in good mutual agreement for log(age Gyr -1 ) > 0.4, 
which indicates that the method based on leaving reddening as 
a free parameter may lead to a partial degeneracy in the result- 
ing ages, since addition of NIR photometry could not break the 
degeneracy completely. As for the equivalent metallicity compar- 
ison, most of our results agree with one another, irrespective of 
the fitting method used, except that some metallicities derived 
from the free-reddening method are much lower than those de- 
rived based on adoption of metallicities determined previously 
in the literature. These results may indeed suffer from the age- 
metallicity-extinction degeneracy. On the other hand, the masses 
derived from both methods are in excellent mutual agreement in 
almost all cases. We also calculated the offsets and la errors of 
the parameters derived from the free- and fixed-reddening fits: 
{E(B - V)frce - E{B - V0 flxcd ) = 0.115 ± 0.233 mag, 
(log Age free -logAge fixcd ) = -0.261±0.628 [yr], ([Fe/H] free - 
[ Fe / H ]flxod> = -0- 133 ± 0- 489 dex ^ <log(M frec /M Q ) - 
log(M flx ed/M )} = -0.039 ± 0.202. The systematic differences 



between the free- and fixed-reddening fits is consistent with zero 
for all parameters. 

Figure [5] shows a similar comparison as Fig. |4] but for the 
GALEV SSP models. In this case, the reddening values derived from 
the free fits, as well as the age estimates, agree better with the litera- 
ture values than in Fig.[4] while the quality and mutual agreement of 
the mass fits is similar as in Fig. [4] However, the metallicity values 
derived from the free-reddening method are systematically higher 
than those based on fits assuming reddening values taken from 
the literature. We thus conclude that the age-metallicity-extinction 
degeneracy affecting the GALEV models is weaker than that for 
the BS-SSP models, but the age-metallicity-extinction degeneracy 
seems worse in this case than that shown in Fig|4] Again, we cal- 
culated the offsets and la uncertainties for the parameters derived 
from the free- and fixed-reddening fits: {E(B — V)i ree — E(B — 
F)fixcd) = 0.028 ± 0.249 mag, (logAge frcc - log Age flxcd ) = 
-0.088±0.740 [yr], ([Fe/H] tree - [Fe/H] fixcd ) = 0.122±0.497 
dex, and (log(M frcc /M ) - log(M flxcd /M )} = 0.032 ± 0.263. 
Similar to Fig. [4] also note that the systematic differences between 
the free- and fixed-reddening fits are consistent with zero for all 
parameters. 

Using Figs [4] and [5] we investigated the extent of the 
age-meta l licity degeneracy for the free-reddening results. As in 
iFanetalJ feoiO). we conclude that there are no significant offsets 
between the fixed- and free-reddening results for any of the ba- 
sic cluster parameters. Therefore, the free-reddening results for the 
310 star clusters with reddening values from the literature will no 
longer be used in the remainder of this paper. 

Figure |6]shows comparisons of the ages, metallicities, redden- 
ing values and masses derived from the BS-SSP and BC03 models. 
The filled data points with error bars in the bottom panels repre- 
sent the mean values for each bin. Each case yields different re- 
sults owing to the different input SSP models and different IMFs 
used. Note that the ages based on the BS-SSP model fits are sys- 
tematically older (by ~0.3 dex) than those derived from the BC03 
models. The results from the BS-SSP models yield higher metallic- 
ity, larger reddening and more massive clusters than those derived 
from the BC03 models. Since the BS-SSP models implicitly as- 
sume that BSs will affect the integrated colours, the clusters appear 
to be younger and more metal poor. Thus, the colours predicted 
by the BS-SSP models are systematically bluer than those result- 
ing from the BC03 model suite for the same age and metallicity. 
Therefore, it is easy to understand that the reddening values de- 
rived from the BS-SSP models are systematically larger than those 
derived from BC03. In addition, the BS-SSP models lead to older 
ages and higher metallicities. Based on Fig. [6] we note that most 
of the outliers are found above the locus of equality, which may 
be caused by relatively large numbers of BSs in those clusters. The 
small number of outliers below the line of equality may be partially 
caused by degeneracies in the fits. Table [Tj presents the systematic 
differences in the ages, metallicities and masses derived from the 
BS-SSP, GALEV and BC03 SSP models. The errors listed in the ta- 
ble are the errors associated with the offsets. They are defined as 
a/yN, where a is the standard deviation and N the number of 
data points. Note that the reddening fits between the BS-SSP and 
BC03 models are consistent with one another, i.e., we do not find 
any systematic differences between the reddening values derived 
from adoption of these models. 

Figure UJis similar to Fig. [6] but uses both the BS-SSP mod- 
els and the GALEV SSP models. The ages resulting from the BS- 
SSP models are systematically older for log(age yr~ x ) > 9.5 (t > 
3.16 Gyr) - see Table []]- while for younger objects the scatter is 
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[Fe/H] free reddenlng '°g(^cl/^©)free reddening 

Figure 4. Comparisons of the reddening values, ages, metallicities and masses derived by adopting reddening values from the literature versus leaving 
extinction as a free parameter for the BS-SSP models. 



Table 1. Mean offsets of the results derived from different pairs of BS-SSP, GALEV and BC03 SSP models for M31 globular-like clusters. The errors quoted 
are the errors associated with the offsets, defined as a/ sfN, where a is the standard deviation and N the number of data points. 



Model suites 


(Alog(Age)) 


<A[Fc/H]> 


(AE(B - V)) 


(Alog(A/ cl )> 




[yr] 


(dex) 


(mag) 


[M ] 


BS-SSP - BC03 


0.307 ± 0.024 


0.050 ± 0.024 


-0.009 ± 0.006 


0.200 ± 0.030 


BS-SSP - GALEV 


0.137 ± 0.022 


0.353 ±0.019 


0.025 ± 0.006 


0.091 ± 0.011 


BC03 - GALEV 


-0.170 ±0.027 


0.303 ± 0.025 


0.035 ± 0.008 


-0.109 ±0.028 



relatively large and no obvious systematic differences are seen. In 
addition, the metallicities resulting from the BS-SSP models are 
~ 0.3 dex higher (more metal-rich) than those derived from the 
GALEV models. This is because the BS-SSP models correct for the 
effects of BSs and predict older ages and higher metallicities than 
the 'standard', single-star models. For the same reason, the result- 
ing reddening values are higher than for the other models (since the 
BS-SSP models predict bluer colours). For the mass estimates, the 
models basically agree with one another. The presence (or absence) 
of BSs will thus significantly affect the fit results. 



Figure[8]is the same as Fig. [6] but the comparisons of the ages, 
metallicities, reddening values and masses are between those de- 
rived from the BC03 models and the GALEV SSP models. The ages 
derived from the BC03 models are slightly younger (by 0.137 dex 
on average) than those derived from the GALEV models, while the 
metallicities derived from the BC03 model suite are significantly 
higher (by ~0.4 dex) than those based on the GALEV set for the 
entire metallicity range. For the masses and reddening values, the 
results are consistent with one another. 

To check whether the differences between the results from 
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Figure 5. As Fig. 4, but for the GALEV SSP models. 



the various models are caused by intrinsic differences among the 
models or owing to poorly understood degeneracies, we fitted the 
metallicities using the BS-SSP and GALEV models, adopting the 
same ages and reddening values as derived from the BC03 mod- 
els. If the metallicities derived from the three models are the same, 
the differences between the results are likely caused by degenera- 
cies in the fits; if not, differences among the models may account 
for the differences in the results. Figure [9] shows comparisons be- 
tween the metallicities derived from the different models but adopt- 
ing the ages and reddening values taken from BC03 models. Note 
that the metallicities are indeed different for the different models 
although the same ages and reddening values have been adopted. 
This thus shows that the differences among the results are likely 
caused by differences among the models. Note that this does not 
necessarily imply that the model differences are not caused by the 
age-metallicity-extinction degeneracy. We calculated the offsets 
between the metallicities derived from the BC03 and BS-SSP mod- 
els for the ages and reddening values adopted based on application 
of the BC03 models: ([Fe/H] BC03 - [Fe/H] BS ) = -0.022±0.269 
dex. For the offsets between the metallicities from the BC03 and 
GALEV models with the same ages and reddening values as for 



BC03 we find ([Fe/H] B C03 - [Fe/H] Ga i cv ) = 0.263 ±0.337 dex. 
Clearly, neither offset is statistically significant. 

Many recent studies show that there are systematic differences 
in the colo urs predicted by diffe rent stellar p opulation synthe- 
sis models. IConrov & GuniJ d2010h note that the lMarastonl 0005!) 
models predict much redde r (J — K), (V — K ) and ( g — r) colours 
for t < 2 Gyr than the iBruzual & Charlotl d2003l) i mode l s, and 
that their own models are more similar to the Maraston d2005l) 
models since the thermally pulsing asymp totic giant branch sec- 
tions of the isochrone sets are comparable. IConrov & GunrJ feoiol) 
obtained old ages for the M3 1 and Galactic massive star clusters 
based on t heir own models, as well as on the IBruzual & Charlotl 
d2003l) and lMaraston] d2005l) models. They found that their FSPS 
models can fit the NIR photometry of star clusters in the Mag- 
ellanic Clouds better than the latest Padova and M05 models. 
The BC03 models also fit these observatio ns quite well, al t hough 
they predict NIR colours that are too blue. ICoMov&Gunnld2010h 
also concluded that all colours - except (J — K) - agree, al- 
though we note that t he scatter in our M31 st a r clust er results is 
fairly large. Recently, iMaraston & Strombackl d201ll) found that 
the MILES-based models predict lower fluxes in the NIR regime 
compared to the STELIB, Pickles or MARCS-based models. In 
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Figure 6. Comparisons of ages, metallicities, reddening values and masses derived from the BS-SSP and BC03 models. In the bottom panels, the filled data 
points with lcr error bars represent the mean values of the binned data. 




Figure 7. As Fig. 6, but for the BS-SSP and GALEV SSP models. 



addition, the MILES-b ased models of IVazdekis et .ail fcOld) and 
IConrov & Gunnl ll2010h agree well with iMaras ton (2009). These 
authors found good agreement of bluer predicted colours, both in 
luminous red galaxies and globular clusters. They also note the 
good agreement between their results from full spec t ral fitt ing and 
colour-magnitud e diagram analys is. [Peacock et al. (201 1) found 
that the iBruzual & Chariot! d2003l) . IMaraston! d2005l) and Padova 
models can fit the M3 1 cluster s well, except the (g — r) colours. 
Maraston & Strombackl J201lh 's MILES models as well as the 
Vazd ekis et alj feoioh models can fit the (g — r) colours of the 
M31 star clusters very well for all metallicities. 



4.2 Age, metallicity and mass distributions 

The cluster age distribution is interesting, because it offers a clue 
to the galaxy's formation history. Figure [10] shows the age distri- 
bution of our sample of globular-like clusters in M31 (bin size: 0.1 
dex). In the top panels, the ages have been derived assuming the 
WMAP age of 13.7 ±0.2 Gyr as upper limit, while in the bottom 
panels the ages have been fitted using the original models without 
imposing an upper age limit, for comparison. Clearly, for the same 
model, the distributions in the top and bottom panels are essentially 
the same for ages < 10 Gyr. However, it seems that there are more 
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Figure 8. As Fig. 6, but for the BC03 and GALEV models. 



star clusters older than the WMAP age in the bottom panels, which 
may be caused by either reddening values that have been underesti- 
mated (i.e., the age-extinction degeneracy) or observational errors. 
It appears that there is a larger fraction of old globular-like clusters 
in the distribution resulting from application of the BS-SSP mod- 
els than in those of the other models, which again confirms that the 
BS-SSP models tend to predict older ages. If we consider the clus- 
ters older than 2 Gyr as the 'old sample', the relative fractions of 
old clusters are 67, 50 and 59 % for the BS-SSP, BC03 and GALEV 
models, respectively. The mean values of the histograms are given 
in Table[2] Note that if we adopt the BC03 models, a large fraction 
of clusters are aged between 1 and 10 Gyr (please see Fig.UO). The 
oldest age estimates are derived from the BS-SSP models, while 
we note that changing the IMF (Kroupa vs. Chabrier) does not ev- 
idently affect the results. There is a large fraction (~ i) of young 
star clus ters (< 2 Gyr) in our distributions. This is similar to the 
result of IWang et aO hold) . 



Table 2. Mean values of the best-fitting age, metallicity and mass distribu- 
tions based on different models. 



The newly derived age distribution of our sample clusters, 
combined with that of the younger clusters reported in previ- 
ous s t udies (e.g-ICaldwefl et alj|2009l ; IWang et alj|2oTol ; iFan et all 
l20ld: IPeacock et al.ll201oF 7 shows evidence of active star forma- 
tion in M31 during the past 2 Gyr. This implies that there may 
have been several star-forming episodes in this period, possibly 
triggered by a major or several minor mergers with other, smaller 
galaxies. The age distribution of the M31 globular-like clusters is 
quite different from that of the Milky Way's GCs, most of which are 
older than 10 Gyr (e.g., IVandenBerg et al.lll996t |Pe Angeli et ajj 
2005). In fact, this might be owing to our sample selection. Our re- 
sults are based on observa tions in the M31 disc, where most of the 
young clusters are located. Illammer et al. I d2007b suggested that the 
Milky Way has had an exceptionally quiet formation history during 
the last 10 Gyr and M31 might have undergone a recent, active 
merger, which may account for the observation that all Galactic 
GCs are old, whereas we find a large fraction of YMCs in M31. 
iMcConnachie et al. (2009) suggested that an encounter between 
M33 and M31 took place a few Gyr ago. The clusters with ages 



Parameter 


BS-SSP 


BC03 


GALEV 


log(age yr" 1 ) 


9.50 


9.19 


9.36 


[Fe/H] (dex) 


-0.475 


-0.525 


-0.828 


log(M cl /M Q ) 


5.16 


4.95 


5.07 



in excess of 10 Gyr in Fig.QJJ]were most likely created during the 
epoch when the galaxy formed, while the young globular-like clus- 
ters might have been created in a number of mergers during the last 
few Gyr or by the postulated recent galactic encounter with M33. 

Figure [TT] shows the resulting metallicity distributions based 
on the different SSP models for a bin size of 0.1 dex. Although 
the upper and lower limits of the metallicity distributions vary for 
the different models, the distributions are different overall and they 
do not exhibit any significant peaks. Table [2] lists the mean val- 
ues of the metallicity distributions. We find that the choice of IMF 
does not significantly affect the results; the average value resulting 
from adoption of the GALEV models represents the lowest metallic- 
ity while that of the BS-SSP models yields the highest metallicity. 
The GALEV models result in a much higher fraction of metal-poor 
([Fe/H] < —1.5) populations compared to both the BS-SSP and 
BC03 models. Since the lower limit in metallicity is much lower 
in the BC03 model suite, the corresponding distribution extends to 
lower metallicities. 

A significant body of recent work focuses on the mass distri- 
bution of the M31 star clusters ( for recent publications see, e.g., 
IFan et al. l l2010l : IWang et alj|2010h . The M/L values can be ob- 
tained from the models as a function of age and metallicity. We cal- 
culated the M/L v values using both the BS-SSP and GALEV SSP 
models, luminosities based on conversion of the V-band fluxes, 
and an M31 distance mo dulus of (m — M)o ~ 24.47 mag 
(McConnac hie et alj2005l) . Figure [T2"lshows the mass distributions 
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Figure 9. Comparisons between the metallicities derived from the different models but for the same ages and reddening values as derived based on the BC03 
models. 



of the M31 star clusters in our sample derived using the BS-SSP, 
GALEV and BC03 models (bin size: 0.2 dex). For each model, the 
mass ranges from ~ 3.4 to ~ 7 dex in units of log(Mo). In ad- 
dition, for comparison we also plot the mass dis tributions of the 
Galactic GCs based on the data of lHarrisI d 19961) (2010 edition; 
see http://www.physics.mcmaster.ca/~harris/mwgc.dati, assuming 
an age of 13 Gyr (but we remind the reader that M/L values are in- 
sensitive to age differences for such old ages). The models applied 
are the same as those used for the M3 1 star clusters. We find that the 
mass distributions of the M3 1 star clusters and the Galactic GCs are 
very similar, irrespective of IMF choice. Note that all mass distri- 
butions have similar peak values, from log(Moi/M©) = 5.0 to 5.5, 
which agrees very well with the universal GC mass function, i.e., 
a Gaussian function with a mean of log(M c i/Mo) ~ 5.2-5.3 and 
a la standard deviation of flWM cl /M ) — 0.5-0.6, as suggested 
bv lParmentier &Gilmorel f2007). Tabled lists the mean values of 
the distributions in Fig.[T2] We find that the BS-SSP models result 
in the most massive mean value, while the BC03 models yield the 
lowest mean mass. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we aimed at exploring the cause(s) of systematic dif- 
ferences among model results based on application of different SSP 
model suites, we have rederived the ages, metallicities, reddening 
values and masses of 445 M3 1 star clusters using \ 2 minimisation 
and adoption of the BS-SSP and GALEV SSP models. Our clus- 
ter sample as well as th e optical broad-b and and NIR photome- 
try used were taken from lFan et al. We also compared the 
SED-matching results for the same sample based on the 'standard' 
single-star BC03 SSP models. 

We compared the U BVRIJHK colours as a function of age 
predicted by different SSP models, specifically the BS-SSP mod- 
els with a Chabrier-like IMF, the BC03 models with a Chabrier 
IMF, and the GALEV models with a Kroupa IMF. We found that the 
BS-corrected models exhibit offsets towards bluer colours for all 
colours, although these offsets tend to become negligible for ages 
< 2 Gyr. However, note that the GALEV models show much red- 
der colours in (V — K) and (J — H) for ages < 2 Gyr, which 
may affect age estimates for young stellar populations. The GALEV 
and BC03 models predict similar colours in (17 — B), (B — V), 
(V — R) and (V — J), while both models exhibit significant offsets 
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Figure 10. Age distributions for different models with different IMFs. Top: Ages derived from the models assuming the WMAPupper age limit. Bottom: Ages 
derived from the original models without imposing an upper age limit. 
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Figure 11. Metallicity distributions for different models. 



in (V — K) and (J — H), especially for young ages. The BS-SSP 
models predict the bluest colours while the GALEV models predict 
the reddest colours. We also found that the GALEV models predict 
the highest M/Ls, while the BS-SSP models predict the lowest 
values. In addition, the differences become larger for older ages. 

Our comparisons show that a reasonable choice of IMF does 
not affect the results significantly and that the main differences 



in the results are caused by intrinsic differences among the SSP 
models. The ages derived from the BC03 models are the youngest 
while those resulting from the BS-SSP models are the oldest. For 
the metallicity, the results of GALEV models are the lowest while 
those from the BS-SSP models are the highest. Finally, the masses 
based on the BC03 models are the lowest and those from the BS- 
SSP models are the highest. 
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Figure 12. Cluster mass distributions for different models. Top: M31 star clusters. Bottom: Galactic star clusters. 



We used photometry in eight filters (UBVRIJHK) for most 
of our sample clusters, where the U band and the NIR JHK filters 
are essential to obtain robust results by means of SED fits (see Sect. 
13 . lb - We also make a case for proper inclusion of the effects of BSs 
in SSP codes. 



ACKNOWLEDGMENTS 

We are indebted to the referee for his/her thoughtful comments 
and insightful suggestions that improved this paper greatly. This re- 
search was supported by the National Natural Science Foundation 
of China (NSFC) under grants No. 11003021, 11043006, 11073001 
and 11073032. 



REFERENCES 

Ahumada, J., Lapasset, E„ 1995, A&As, 109, 375 
Anders P., Fritze-v. Alvensleben U., 2003, A&A, 401, 1063 
Anders, P., Fritze-v. Alvensleben, U., de Grijs, R., 2003a, ASPC, 
296, 567 

Anders, P., Fritze-v. Alvensleben, U., 2003b, A&A, 401, 1063 
Anders P., de Grijs R., Fritze-v. Alvensleben U., & Bissantz N, 

2004a, MNRAS, 347, 17 
Anders, P., Bissantz, N, Fritze-v. Alvensleben, U. & de Grijs R., 

2004b, MNRAS, 347, 196 



Arimoto, N., 1996, From Stars to Galaxies, Leitherer C, Fritze- 
v. Alvensleben U., Huchra J., eds., ASP Conf. Ser., (ASP: San 
Francisco), 98, p287 

Ascenso, J., Alves, J., Beletsky, Y., & Lago, M. T. V. T. 2007a, 
A&A, 466, 137 

Ascenso, J., Alves, J., Vicente, S., & Lago, M. T. V. T. 2007b, 

A&A, 476, 199 
Barmby, P., et al. 2000, AJ, 1 19, 727 
Barmby, P., & Huchra, J. P. 2001, AJ, 122, 2458 
Barmby, P., et al. 2009, AJ, 138, 1667 

Bastian, N., Gieles, M., Lamers, H. J. G. L. M., Scheepmaker, R. 

A., de Grijs, R. 2005a, A&A, 431, 905 
Bastian, N., Gieles, M., Efremov, Yu. N., Lamers, H. J. G. L. M. 

2005b, A&A, 443, 79 
Battistini, P., Bonoli, F, Braccesi, A., Fusi-Pecci, F, Malagnini, 

M. L., & Marano, B. 1980, A&As, 42, 357 
Battistini, P., Bonoli, F, Braccesi, A., Federici, L., Fusi Pecci, F, 

Marano, B., & Borngen, F. 1987, A&As, 67, 447 
Battistini, P., Bonoli, F, Casavecchia, M., Ciotti, L., Federici, L., 

& Fusi Pecci, F. 1993, A&A, 272, 77 
Brazual, G, & Chariot, S. 1993, ApJ, 405, 538 
Bruzual, G, & Chariot, S. 2003, MNRAS, 344, 1000 (BC03) 
Caldwell, N., Harding, P., Morrison, H., Rose, J. A., Schiavon, R., 

& Kriessler, J. 2009, AJ, 137, 94 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 

Carpenter, J. M., Hillenbrand, L. A., & Skrutskie, M. F. 2001, AJ, 



Systematic differences in SSP models: M31 star clusters 15 



121, 3160 

Cenarro, A. J., Cervantes, J. L., Beasley, M. A., Marin-Franch, A., 

& Vazdekis, A., 2008, ApJ, 689, L29 
Chabrier, G. 2003, PASP, 115, 763 
Chariot, S., & Bruzual, G. 1991, ApJ, 367, 126 
Conroy, C, Gunn, J. E., White, M. 2009, ApJ, 699, 486 
Conroy, C, & Gunn, J. E. 2010, ApJ, 712, 833 
Crampton, D., Cowley, A. P., Schade, D., & Chayer, P. 1985, ApJ, 

288, 494 

De Angeli, E, et al. 2005, AJ, 130, 1 16 

de Grijs, R., Bastian, N., & Lamers, H. J. G. L. M. 2003, MNRAS, 
340, 197 

de Grijs, R., Fritze-v. Alvensleben, U., Anders, P., Gallagher, J. S., 
Bastian, N., Taylor, V. A., Windhorst, R. A. 2003, MNRAS, 342, 
259 

de Grijs, R., et al. 2004, MNRAS, 352, 263 

de Grijs, R., et al. 2005, MNRAS, 359, 874 

de Grijs, R., & Anders, P. 2006, MNRAS, 366, 295 

Fan, X., et al. 1996, AJ, 112, 628 

Fan, Z., Ma, J., de Grijs, R., Yang, Y, & Zhou, X. 2006, MNRAS, 
371, 1648 

Fan, Z., Ma, J., de Grijs, R., & Zhou, X. 2008, MNRAS, 385, 
1973 

Fan, Z., de Grijs, R., & Zhou, X. 2010, ApJ, 725, 200 
Fioc, M„ & Rocca-Volmerange, B. 1997, A&A, 326, 950 
Frayn, C. M., & Gilmore, G. F. 2002, MNRAS, 337, 445 
Fritze-v. Alvensleben, U., Bicker, J. 2006, A&A, 454, 67 
Galleti, S., Federici, L., Bellazzini, M., Fusi Pecci, E, & Macrina, 

S. 2004, A&A, 416, 917 
Galleti, S., Federici, L., Bellazzini, M., Buzzoni, A., & Fusi Pecci, 

F. 2006, A&A, 456, 985 
Galleti, S., Bellazzini, M., Federici, L., Buzzoni, A., & Fusi Pecci, 

F. 2007, A&A, 471, 127 
Galleti, S., Bellazzini, M., Buzzoni, A., Federici, L., Fusi Pecci, 

F. 2009, A&A, 508, 1285 
Gieles, M., Bastian, N., Lamers, H. J. G. L. M., Mout, J. N., 2005, 

A&A, 441, 949 

Hammer, F, Puech, M., Chemin, L., Flores, H., & Lehnert, M. D. 

2007, ApJ, 662, 322 
Harris, W. E. 1996, AJ, 112, 1487. 

Hodge, P., Krienke, O. K., Bianchi, L., Massey, P., Olsen, K. 

2010, PASP, 122, 745 
Hubble, E. 1932, ApJ, 76, 44 

Hunter, D. A., O'Connell, R. W., Gallagher, J. S., Ill 1994, AJ, 
108, 84 

Jiang, L., Ma, J., Zhou, X., Chen, J., Wu, H., & Jiang, Z. 2003, 
AJ, 125, 727 

Kaviraj, S., Rey, S. C, Rich, R. M., Lee, Y. W., Yoon, S. J., Yi, 
S. K. 2007, MNRAS, 381, 74 

Kotulla, R, Fritze, U., Weilbacher, P., & Anders, P. 2009, MN- 
RAS, 396, 462. 

Krienke, O. K., & Hodge, P. W. 2007, PASP, 1 19, 7 

Kroupa, P. 2001, MNRAS, 322, 231 

Kurucz, R. L. 1992, in IAU Symp. 149, The Stellar Populations of 
Galaxies, ed. B. Barbuy & A. Renzini (Dordrecht: Kluwer), 225 
First citation in article 

Lejeune, T, Cuisinier, E, & Buser, R. 1997, A&As, 125, 229 

Lejeune, T, Cuisinier, E, & Buser, R. 1998, A&As, 130, 65 

Leitherer, C, & Heckman, T. M. 1995, ApJS, 96, 9 

Ma, J., et al. 2007, ApJ, 659, 359 

Ma, J., et al. 2009, AJ, 137, 4884 

Macri L. M. 2001, ApJ, 549, 721 



Mafz Apellaniz J., 2009, ApSS, 324, 95 
Maraston, C. 1998, MNRAS, 300, 872 
Maraston, C. 2005, MNRAS, 362, 799 

Maraston, C, Nieves Colmenarez, L., Bender, R., Thomas, D. 

2009, A&A, 493, 425 
Maraston, C, & Stromback, G. 2011, MNRAS, 418, 2785 
Massey, P., Olsen, K. A. G, Hodge, P. W., Strong, S. B., Jacoby, 

G. H., Schlingman, W., & Smith, R. C. 2006, AJ, 131, 2478 
McConnachie, A. W., et al. 2005, MNRAS, 356, 979 
McConnachie, A. W., et al. 2009, Nature, 461, 66 
McCrady, N. 2009, Ap&SS, 324, 109 

Moll, S. L., de Grijs, R., Mengel, S., & Smith, L. J. 2008, ASPC, 
388,411 

O'Connell, R. W., Mangano, J. J. 1978, ApJ, 221, 62 
Parmentier, G., & Gilmore, G. 2007, MNRAS, 377, 352 
Peacock, M. B., Maccarone, T. J., Knigge, C, Kundu, A., Waters, 

C. Z., Zepf, S. E., & Zurek, D. R. 2010, MNRAS, 402, 803 
Peacock, M. B., Zepf, S. E., Maccarone, T. J., & Kundu, A. 2011, 

MNRAS, 737, 5 
Perina, S., et al. 2010, A&A, 494, 933 
Perina, S., et al. 2010, A&A, 511, 23 
Piotto, G, et al. 2002, A&A, 391, 945 
SalpeterE. E. 1955, ApJ, 121, 161 

Sargent, W. L. W., Kowal, C. T, Hartwick, F. D. A., & van den 

Bergh, S. 1977, AJ, 82, 947 
Schulz, J., Fritze-v. Alvensleben, U., Moller, C. S., Fricke, K. J. 

2002, A&A, 392, 1 
Schulz, J., Fritze-v. Alvensleben, U., Fricke, K. J. 2003, A&A, 

398, 89 

Smith, L. J., Kowal, et al. 2007, ApJ, 667, L145 
Spergel, D. N, et al. 2003, ApJS, 148, 175 
Stanek, K. Z., & Garnavich, P. M. 1998, ApJ, 503, 131 
VandenBerg, D. A., Bolte, M., Stetson, P. B. 1996, ARA&A, 34, 
461 

Vazdekis, A. 1999, ApJ, 513, 224 

Vazdekis, A., et al. 2010, MNRAS, 404, 1639 

Vetesnik, M. 1962, Bulletin of the Astronomical Institutes of 

Czechoslovakia, 13, 180 
Wang, S., Fan, Z., Ma, J., de Grijs, R., & Zhou, X. 2010, AJ, 139, 

1438 

Whitmore, B., et al. 2005, AJ, 130, 2104 
Worthey, G. 1994, ApJs, 95, 107 

Wu, H., Shao, Z. Y, Mo, H. J., Xia, X. Y, & Deng, Z. G. 2005, 

ApJ, 622, 244 
Xin, Y, & Deng, L., 2005, ApJ, 619, 824 
Xin, Y, Deng, L., & Han, Z. W„ 2007, ApJ, 660, 319 
Xin, Y, Deng, L. C, de Grijs, R., & Kroupa, P. 2011, MNRAS, 

411,761 

Yi, S. K, Kim, Y. C, & Demarque, P. 2003, ApJS, 144, 259 



